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A method (Fig. 1) for creating a frequency or time 
domain semblance for use in conjunction with acoustic logging 
tools (100, 130, 140, 145). The frequency domain semblance 
(Fig. 2F) may be obtained by transforming an acoustic signal 
(150) received at multiple depths into the frequency domain, 
combining the received waveforms corresponding to the different 
depths, and expressing the result in a graph with slowness and 
frequency axes. Another aspect of the invention is the treatment 
of two or more time domain semblances (Fig. 2B) as probability 
density functions of the slowness for an acoustic signal. This 
enables the combination of time domain semblances from the 
same depth in the wellbore. Once combined the time domain 
semblances more accurately depict the slowness of an acoustic 
wave through the formation at the selected depth. In addition, 
a related self-adaptation method to compress the waveform data 
downhole for storage or transmission is also disclosed. 
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PROCESSING FOR SONIC WAVEFORMS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application relates to provisional application serial number 60/128,912 filed April 12, 

1999. 

5 STATEMENT REGARDING FEDERALLY SPONSORED 

RESEARCH OR DEVELOPMENT 

Not Applicable. 

BACKGROUND OF THE INVENTION 
To determine whether a particular rock formation contains oil or other hydrocarbons, it 

1 0 must have certain properties. Acoustic tools, also known as sonic tools, are one way to investigate 
a rock formation around a weilbore. 

As shown in Figure 1, an acoustic tool 100 may be part of a bottomhole assembly in a drill 
string 110 drilling through a rock formation 120. Alternately, such an acoustic tool may be part of 
a wireline device. The acoustic tool may include an acoustic transmitter 130 and a set of acoustic 

15 receivers 140, 145. When the acoustic tool is placed in the weilbore drilled through a rock 
formation, the tool 100 transmits a signal 150 from its transmitter 130. This signal travels through 
the rock formation 120 and arrives at the receivers 140, 145. The signal is detected at the set of 
receivers as a series of waveforms 155 as generally shown in Figure 2A. One parameter of 
particular interest is an acoustic signal's speed through the surrounding formation. Speed can also 

20 be expressed as slowness, the inverse of speed Such information is then used to infer whether 
hydrocarbons are present in the rock formation around the weilbore. As shown in Figure 2A, a set 
of seven receivers generates seven different detected waveforms or "channels." Based on the time 
delay of the waveforms or channels, the slowness of the signal through the rock formation and 
other characteristics of the rock formation can be determined. 

25 The data from acoustic tools may be placed on a graph for simplified interpretation. One 

known method to analyze the waveforms is a time domain semblance as shown in Figure 2B. The 
time domain semblance for a set of waveforms is obtained by stacking the waveforms or channels 
through a range of slowness values. A range of semblance values are then assigned to the stacked 
waveforms, with a higher semblance value corresponding to a higher degree of "fit" among all the 

30 waveforms. 

This resulting time domain semblance is often a color-coded "map" drawn on a time axis 
and a slowness axis. The locations of "peaks" on die map, shown by high semblance values, 
indicate the estimate of the slowness (the inverse of velocity) of the received signal from the 
transmitter and the time of arrival for the signal at the last receiver (although any receiver could be 
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chosen to correlate with the time axis). For example, a low semblance value may be represented 
by a blue color, a medium semblance value may be represented by a yellow color, and a high 
semblance value may be represented by a red color. 

Complicating this analysis is that after being generated by the acoustic transmitter, and 
5 depending upon the frequency of the acoustic signal and the characteristics of a rock formation, an 
acoustic signal from the transmitter may excite a variety of types of secondary acoustic waves. 
These types of secondary acoustic waves include compressional waves, shear waves, and Stonely 
waves. Each of these may then be indicated by a different peak on the time domain semblance. In 
addition to showing the location of received waves on the semblance, the time domain semblance 

10 gives an indication of the intensity of the received wave and a "shape" corresponding to each 
received wave. Multiple time domain semblances are then be used to create a log or other 
correlation between depth and slowness. 

However, semblance graphs have numerous drawbacks, including the presence of aliases, 
also known as shadows, on the slowness/time semblance graph. These aliases or shadows are 

15 spots or peaks on the semblance graph that indicate that a wave was received at the acoustic 
receivers, but in reality these peaks are phantom "shadows" that exist simply because of problems 
in the derivation of the semblance graph. This creates uncertainty regarding the analysis of the 
traveling wave, and can therefore lead to errors in the analysis of the formation. 

In addition, a "smeared" or spread out peak on the semblance graph can also arise from 

20 unusually high dispersion of the acoustic signal in a rock formation, such as happens in a very soft 
rock formation. Similar smearing or diffusion of a peak can result from problems with the acoustic 
tool. Thus, the time domain semblance has difficulty distinguishing between an improperly 
operating tool and a rock formation that disperses acoustic waves to an unusually high degree. 

Further, smearing and other spread-out peaks are a problem because of the difficulty of 

25 finding the "true" location of a wave's slowness. Thus, the smearing of a peak decreases the 
accuracy of the measurements and increases the chances that a mistake will be made. A solution is 
needed that reduces or solves these problems. 

Another problem in acoustic logging is the enormous amount of data that must be collected 
while the tool is downhole. Memory carried within the acoustic tool is expensive and should not 

30 be wasted storing unnecessary data. Thus, downhole compression of the waveform data is 
necessary. Nonetheless, even more efficient compression of the data is desired to further conserve 
downhole memory. 

Yet another problem in logging is to compute the slowness of different types of waves at 
various depths in the wellbore. One attempted solution to this problem is to store waveform data 
35 downhole and analyze the data at the surface. However, the analysis of the data at the surface to 
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find peaks in semblance and other characteristics of the received waves is very slow and wasted 
time in the field is expensive. This solution is therefore undesirable. Another approach to this 
problem is to program a downhole processor to identify waves and their slowness downhole. 
However, while this approach substantially speeds the analysis (by utilizing "dead time" between 
5 consecutive acquisitions) it is inaccurate because of the complexities in analyzing wellbore data. 
Thus, a better solution is required. Ideally, this solution would combine the efficiency of 
performing some processing downhole with the accuracy of the surface analysis. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a side view of a prior art acoustic tool in a wellbore; 
10 Figure 2 A is a set of waveform time series showing the reception of a traveling waveform 

at a spaced set of receivers; 

Figure 2B is a time domain semblance of the waveforms of Figure 2 A; 
Figure 2C is an outline view of Figure 2B 

Figure 2D is the frequency response for the seven channels shown in Figure 2A; 
1 5 Figure 2E is a frequency semblance corresponding to the time domain semblance of Figure 

2B; 

Figure 2F is an outline view of Figure 2E; 

Figure 2G is a frequency semblance corresponding to that shown in Figure 2E but over a 
broader range of frequencies; 
20 Figure 3 is a novel acoustic tool including more than one set of acoustic receivers spaced 

around the circumference of the tool; 

Figure 4A is a first frequency domain semblance; 
Figure 4B is an outline view of Figure 4 A; 
Figure 4C is a second frequency domain semblance; 
25 Figure 4D is an outline view of Figure 4C; 

Figure 4E is a third frequency domain semblance; 

Figure 5 is an analog waveform illustrating a self-adaptable compression technique; 
Figure 6A is a time semblance for a set of data; 
Figure 6B is an outline view of Figure 6A; 
30 Figure 6C is a frequency semblance for the set of data of Figure 6A; 

Figure 6D is an outline view of Figure 6C; 

Figure 6E is the frequency semblance of Figure 6C over a broader range of slowness 

values; 

Figure 6F is an outline view of Figure 6E. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
One aspect of the invention is a novel approach to interpret time domain semblance graphs. 
In particular, the invention interprets semblance graphs as probability density functions (layouts). 
This allows a series of semblance graphs to be combined to yield a more accurate estimate of 
5 waveform slowness and other characteristics of the rock formation around the wellbore. 

Inaccuracy and smearing of a peak on a semblance graph may occur because an acoustic 
tool is not centered in the wellbore. The invention solves this problem first by a novel sonic tool 
design. Referring to Figure 3, it can be seen that a sonic tool according to this aspect of the 
invention includes more than one set of acoustic receivers spaced around the sides of the acoustic 

10 tool Multiple acoustic transmitters may also be spaced evenly around the outer perimeter of the 
tool. The two or more transmitters would then be fired at the same time. This results in more than 
one semblance graph for the same depth. Alternately, acoustic data for the same depth could be 
obtained by other methods or devices. 

By interpreting semblance graphs for the same depth as not only broad indications of time 

15 and slowness, but instead as probability density functions, two or more graphs may be combined 
more effectively to result in a more accurate estimate of the wave slowness and time at a particular 
depth. For example, one semblance graph from a particular depth may have a more spread-out 
peak for a particular slowness and time, while the additional semblance graph from that depth may 
have a more concentrated peak corresponding to that formation. To estimate most closely the 

20 actual slowness and time for the formation, the two graphs should be treated as a density function 
that weights the peaks as probabilities. The semblance surfaces (Le. array of semblance values) 
corresponding to the same depth are then multiplied together to form a maximum likelihood 
estimate of the slowness of the formation. For example, the semblance value at a particular 
slowness and time from one acoustic signal will be multiplied with the semblance value at the 

25 same slowness and time from a second acoustic signal (and a third, etc). The slowness 
corresponding to the peak value of the maximum likelihood estimate is used as the estimated 
slowness for the formation. Thus, if one of the measurements has a higher semblance value for a 
certain slowness, this slowness value results in being given a greater weight in the final 
combination of the values. A diffuse or spread out peak is considered less well determined, will 

30 commonly have a lower semblance value and ultimately has less impact or weighting in the final 
combination. Other forms of maximum likelihood estimates (other than straight multiplication) 
can be utilized as well. Of course, it would be understood by one of ordinary skill upon reading 
this disclosure that calculations may be performed directly on the anay of semblance data without 
any need to first construct a color-coded graph. 
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Another aspect of the invention is the creation of a frequency domain semblance. Just as a 
time domain semblance is a time-slowness graph portraying a waveform's slowness through a 
particular formation, the frequency domain sequence is a frequency-slowness graph that depicts 
that waveform's slowness. Like the time domain semblance, a frequency domain semblance may 
5 be color coded, with blue colors represented low semblance values, yellow colors representing 
medium semblance values, and red colors representing high semblance values. Such a frequency 
domain semblance can be used either on its own to determine the characteristics of a rock 
formation, or can be used in conjunction with a time domain semblance to increase the accuracy 
and reliability of measurements regarding the rock formation. 
10 Referring generally to Figure 2A, a series of waveforms arrive at a set of receivers, for 

example, seven. As can be seen, the acoustic signal arrives initially at the closest receiver (channel 
1 for the example) and last at the receiver furthest away from the transmitter (channel 7 for the 
example). 

A typical tiYne domain semblance 200 for these waveforms is shown in Figure 2B. Figure 
15 2C is a traced outline of the semblance shown in Figure 2B, including labeled peaks 210, 220, 230, 
240. In both Figures 2B and 2C, a peak 210 of relatively high intensity or semblance value exists 
somewhere around 100 microseconds/ft, with another peak 220 at about 225 microseconds/ft. It 
can be seen by reference to Figures 2B and 2C that some uncertainty exists on the semblance graph 
as to the arrival time and speed of the received signal. In addition, a question exists as to how to 
20 interpret what type of wave each peak represents and whether a particular peak is simply a shadow 
or alias of some other peak. Until now, the best way to deteraiine what each peak represents was 
through either an experienced eye or by use of a peak-picking software program, which may be 
subject to error. 

Referring to Figure 2E 5 the frequency semblance corresponding to Figure 2B is shown. 

25 Figure 2F is an outline trace of Figure 2E. Referring to Figures 2E and 2F, there are three notable 
peaks of high semblance value. A first peak 260 is at about 100 microseconds per foot at a 
frequency of somewhere between 1 5 and 20 kilohertz. A second peak is present at about the same 
frequency as the first peak, but has a slowness of about 225 microseconds per foot. This second 
peak 270 is at about the same frequency as the first peak 260, but is distinctly angled from 

30 horizontal. A third peak 280 also exists at about 225 microseconds per foot but at a lower 
frequency somewhere between 5 and 10 kilohertz. 

Figures 2E and 2F can be used to identify an alias or shadow in Figure 2B. In particular, 
the second peak 270 in Figure 2D is a shadow of the first peak 260. This is apparent because of the 
slope or angle of the second peak 270 as compared to the first peak 260 that has little or no slope. 

35 To assist in determining exactly which peak this corresponds to in the time semblance of Figure 
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2B, for example a band stop filter may be applied to the time domain semblance to filter out the 
shadow peak on the time domain semblance. Thus, peak 220 on Figures 2B and 2C can be 
identified as a shadow or alias. 

The angled slope of aliases in the frequency semblance has additional advantages. Figure 
5 2G is a frequency semblance corresponding to about the same depth as Figure 2B, but displays a 
broader range of frequencies. Looking at the broader range of frequencies of Figure 2G, it can be 
seen that the slope of the aliases results in a number of angled peaks that point at the true or actual 
peak. 

In addition, the third peak 280 on Figures 2E and 2F is of interest. The frequency of a peak 
10 tends to indicate how the wave propagated, and thus what type of wave it is. For example, 
compressional waves tend to have a higher frequency, shear waves tend to have lower frequencies, 
and Stonely waves tend to have even lower frequencies. Therefore, the low frequency of peak 280 
tends to indicate that it is not a compressional wave. Moreover, the shape of the third peak on 
Figure 2D indicates that it corresponds to a shear wave. In particular, shear wave peaks tend to 
15 droop or bend down from lower to higher frequencies (left to right in the figure). Other shapes 
may also correspond to a particular type of wave or other characteristic. Consequently, the shape 
of the peaks in a frequency semblance can also be of interest. Of course, although a color-coded 
frequency semblance graph is helpful because it visually and dramatically illustrates the 
characteristics of the received acoustic waveform, these frequency identifications described herein 
20 can also be made directly to the observed data without the need to construct a frequency 
semblance, such as by an automated peak picker that determine the frequency and slowness value 
for maximum semblance, as well as extent and geometric shape of the areas of high semblance 
values. 

Another problem that is solved by a frequency domain semblance is not knowing whether a 
25 diffuse peak is due to dispersion of the acoustic signal in a rock formation or due to tool error. A 
frequency domain semblance resolves this ambiguity because on the frequency domain semblance 
a received waveform in a dispersive formation has a signature shape. In particular, as shown in 
Figures 6A-6F, a peak from a wave that traveled through a dispersive formation can be spread out 
along the frequency axis. Figures 6A and 6B depict a time semblance for a particular set of data, 
30 including a first peak 610 and a second peak 620. Figures 6C and 6D are the corresponding 
frequency semblance including peaks 630 and 640. Peak 630 is the "real" or "actual" peak, 
whereas peak 640 is merely an alias or shadow. 

As can be seen by reference to these Figures, the spread-out shape of the "real" peak 630 in 
the frequency semblance tends upward at higher frequencies when the acoustic signal has 
35 dispersed through a rock formation. This indicates to an operator that the smeared peak present in 
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10 



15 



20 



25 



a time domain sequence is not due to tool error, but instead is due to dispersion of the acoustic 
signal into the rock formation. Figures 6E and 6F are frequency semblances for the same data, but 
over a larger range of slownesses. These Figures include peaks 650, 660, and 670. It can be seen 
that peak 630 at a slowness of around 60-100 microseconds/foot is the "real" peak, and not an 
alias, because of the way peaks 650 and 670 t4 point" toward peak 660, which is also at a slowness 
of around 60-100 microseconds/foot. Thus, peak 660 corresponds to peak 630 and consequently 
both of peaks 660 and 630 are the "real" peak. Using frequency semblances, as shown above, it is 
possible to determine the dispersive nature of a formation. As with time domain sequences, a log 
or other correlation of depth and slowness may be created. This in turn helps to infer the actual 
lithology of the formation. This is another strength that frequency semblance brings to sonic tool 
data processing. 

The frequency domain semblance also has another advantage. As stated above, a 
frequency filter may be applied to the waveforms that make up the time domain semblance to 
remove a shadow or alias. For example, Figures 4A and 4B represent a time domain semblance 
for a particular set of received waveforms. Labeled in Figure 4B is a first peak 410, a second peak 
420, and a third peak 430. Figures 4C and 4D are the corresponding frequency domain semblance. 
Included in Figure 4D are peaks 440 and 450. As can be seen by the slope of peak 450, second 
peak 420 in Figure 4A is a shadow of first peak 4 1 0. When a frequency filter is applied to the time 
domain semblance shown in Figure 4A, it results in the time domain semblance of Figure 4E. As 
can be seen, the peak corresponding to third peak 430 is now stronger and more easily identified 
because of the noise that was removed by use of the frequency filter. In addition, Figure 4E does 
not include the difiuse peak 410 because this shape and slowness of peak has been filtered out as a 
known consequence of resonating the tool body with an acoustic transmitter. 

The data displayed in a frequency domain semblance can be created by the following 
mathematical formula: 



F(x) = Fourier transform; 

fj(t) = time series values for a particular receiver or channel j ; 

d - distance from a constant point (usually the location of the transmitter); 

j = receiver or channel number, 

s = slowness, and 




Where, 
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co = frequency. 

Thus, for any desired frequency and slowness, a frequency semblance value (FS) may be 
derived. 

As can be appreciated by one of ordinary skill in the art upon review of this equation, the 
5 wave forms from the receivers on the sonic logging tool are transformed by use of a Fourier 
transform into the frequency domain to derive fj(t). This results in frequency waveforms as can be 
seen in Figure 2C. The data from each different channel (i.e. receivers), fftt), are multiplied by a 
phase factor e"* djs that corresponds to the different formation slowness (velocity) with which the 
signal propagated through the formation. The phase factors are determined according to the above 
10 equation. 

The signals from each channel are summed for every frequency and slowness that are of 
interest The absolute value is taken of the sum corresponding to one frequency and one slowness. 
This sum is then squared. Alternatively, the squared value may be found by multiplying the 
Fourier sums with its conjugate. Either way, this squared value is then divided by the sum of the 
15 squared absolute value for which the semblance has been computed, as shown in the equation. If 
the phase of the signals from the different channels at a certain frequency coincide for a certain 
slowness phase factor, the value of the frequency semblance will be large and indicate the phase 
velocity of the signal at the analyzed frequency. 

The invention also encompasses a self-adaptable data compression method for acoustic 
20 data. The invention selects an optimum compression technique from a series of available 
compression techniques to minimize the amount of data required to be stored or transmitted. 

Figure 5 shows an analog waveform 500, and a number of sample points 510-518 
corresponding to the digitization of the waveform. For example, when digitized a signal may be 
sampled 10 times per wavelength. As can be seen each consecutive sample point 510-518 has a 
25 step difference 520 from the previous sample point. 

It is not necessary to save the actual value of every sample point 510-518. While the actual 
value of the first point 510 must be saved, the value of the second point 51 1 may be stored as the 
magnitude of step difference 520 from the first point 510. Similarly, the value of the third point 
5 12 may be stored as the magnitude of the difference from the second sample point 511. Thus, the 
30 waveform may be stored as a series of numbers (bits), with the actual value of the first sample 
point being recorded, and followed by the step difference for each consecutive sample point. In 
computer readable form, a bit series of a predetermined number of bits is assigned to record each 
step difference. 

The difference between certain sample points, such as points 5 1 3 and 5 1 4, may be too great 
35 to be recorded by the number of bits assigned to each step difference (ie., the bit series). The 
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magnitude of this difference is labeled 525 in Figure 5, and the maximum difference recordable by 
the predetermined number of bits dedicated to each step difference is labeled as 530. In such an 
instance, an upshift or downshift may be inserted in the data to represent the upshift or downshift. 
For example, for a 4-bit step size an upshift may be defined as 1 1 1 1 and a downshift by 0000. 
5 Consequently, because the upshift and downshift steps occupy two of the 2 n possible combinations, 
a difference amount of (2 n - 2) is the largest difference step 530 that can be represented by the n-bit 
series. 

This 2 n - 2 size difference step can be recorded according to various protocols. For 
example, the middle (or any other defined bit sequence) of the bit series may correspond to a 

1 0 difference step of zero, with a bit series such as 000 1 being a down step and 1110 being an up step. 
Alternatively, a particular bit pattern may indicate a step direction change, with all the consecutive 
step differences being in the same direction (until the "change step direction" bit pattern is once 
again encountered). The size of the difference step in this instance would of course be 2 n - 2. This 
general encoding formula may be referred to as a bit encoding protocol. Other storage protocols 

15 are within the skill of the ordinary artisan. Depending on the particular waveform being 
compressed, and other factors, one or the other of these may be more efficient. 

Another difficulty is determining the ideal number of bits to dedicate for each step 
difference. If too many bits are dedicated to define the size of a step, so that for example the step 
difference between sample points 513 and 514 can be recorded with only one bit series, then often 

20 other step differences along the waveform do not utilize all the bits in the bit series. Thus, bits go 
unused and memory is wasted. If too few bits are dedicated to define the size of a step, then 
insertion of an upshift or downshift in the data sequence will be common and this will result in a 
less efficient than ideal compression. This problem is particularly thorny in the context of acoustic 
wave receivers spaced a range of distances away from an acoustic transmitter because the 

25 waveforms from the more distant receivers tend to have smaller amplitudes, and therefore smaller 
step differences than the waveforms from the closer receivers. Thus, the ideal compression size 
may change depending on the receiver (and waveform) at issue. 

The solution to this dilemma lies in self-adaptation. Before being stored, the amount of 
data needed by each waveform for storage is computed by a micro-processor or the like based on a 

30 range of possible bit values and storage protocols. For example, the amount of memory required to 
store a received waveform may be determined for a 3-bit difference storage scheme, a 4-bit 
difference storage scheme, a 5-bit difference storage scheme, etc. up until a 12-bit storage scheme. 
Of course, the range of bit sizes may be smaller or greater depending on the circumstances. The 
amount of memory required to store a received waveform may also be determined based on the 

35 middle of the bit series corresponding to a bit step of zero, in addition to being computed according 
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to various other methodology, including those described above. The microprocessor then selects 
the step bit size that will ultimately use the fewest number of bits after compression. Afterwards, 
the microprocessor compresses the waveform according to that step bit size. Further, there is no 
need to run actual compression if the microprocessor determines that even the best compression 
5 does not save any memory. 

Those skilled in the art will recognize that many variations on this theme exist. For 
example, compression is not limited simply to the difference between consecutive sample points 
on a waveform. The difference of the differences may also be used, or the difference of the 
differences of the differences, etc. In formula, this may be expressed as: 
10 d(n) = s(n+l)-s(n); 

d(n+l) = s(n+2)-s(n+l); 

dd(n) = d(n+l)-d(n); 

dd(n+l) = d(n+2)-d(n+l); 

ddd(n) = dd(n+l)-dd(n); 
15 Where, 

s(n) = the waveform signal at location n 

d(n) = the difference of consecutive sample points 

dd(n) = the difference of consecutive differences, 

etc. 

20 If the difference of the differences is utilized, then the value of the first difference should be 

stored in addition to the value of the first point in the waveform. This logic then continues if the 
difference of the differences of the differences is used, etc. Some or each of these difference of the 
differences approaches, and other techniques to compress data, can be employed by the downhole 
microprocessor to self-adapt and establish the ideal compression for any particular waveform. 

25 Once the waveform has been compressed, it may be stored or transmitted to the surface. 

After the compression and storage of the data, a downhole processor or processors may 
also start to analyze the data. Because drilling continues for long periods of time, a substantial 
amount of analysis may be completed downhole. This allows the downhole processor to determine 
the characteristics of the waveforms. In particular, the processing can begin the evaluation of 

30 semblance maps for peaks, including their location, the semblance or intensity value, the shape of 
the peak, and the energy of the peak. Thus, a large portion of the time required to establish these 
parameters is eliminated from uphole analysis. On the other hand, the invention does not attempt 
to identify the wave types or other characteristics of the waveforms until the data has been brought 
uphole. This avoids many of the errors associated with previous attempts at downhole waveform 

35 analysis. 
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While preferred embodiments of this invention have been shown and described, 
modifications thereof can be made by one skilled in the art without departing from the spirit or 
teaching of this invention. The embodiments described herein are exemplary only and are not 
limiting. Many variations and modifications of the system and apparatus are possible and are within 
5 the scope of the invention. 
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CLAIMS 

1 . A method to measure formation properties, comprising: 

(a) generating a first acoustic signal at a first depth in a borehole; 

(b) receiving said first acoustic signal at a first set of acoustic receivers 
5 proximate said first depth to generate a first series of received acoustic waveforms; 

(c) calculating from said first series of received acoustic waves a probability 
density function of slowness and time for said first acoustic signal; 

(d) generating a second acoustic signal at said first depth in said borehole; 

(e) receiving said second acoustic signal at a second set of acoustic receivers 
10 proximate said first depth to generate a second series of received acoustic waveforms, 

wherein said second set of acoustic receivers may be different from or the same as said first 
set of acoustic receivers; 

(f) calculating from said second series of received acoustic waves a probability 
density function of slowness and time for said second acoustic signal; 

15 (g) combining said probability density functions of at least said first and second 

series of acoustic waves to derive a most probably measured value formation characteristic 
for fonnation near said borehole at said first depth. 

2. The method of claim 1 , wherein said most probably measured value includes a most 
likely slowness value. 

20 3. The method of claim 1, wherein said calculating steps include stacking a series of 

received acoustic waveforms. 

4. The method of claim 1, wherein said combining step includes multiplying said 

probability density functions. 

5 . A frequency domain semblances method comprising: 

25 (a) transforming data from a first acoustic signal into the frequency domain; 

(b) calculating a frequency semblance value from said data, said frequency 
semblance value corresponding to a first frequency and a first slowness for said data. 

6. The frequency domain method of claim 5, wherein said calculating step is 
conducted to the equation: 

2 

?F(fj(t))e-* 



30 FS(co,s)= 



?|F(f i (t)> 



Where, 

FS(©,s) = frequency semblance value; 
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F(x) = Fourier transform; 

fj(t) = time series values for a particular receiver or channel j; 
d = distance from a constant point (usually the location of the transmitter); 
j = receiver or channel number; 
5 s = slowness, and 

(0 - frequency. 

7. The frequency domain method of claim 5, wherein said transforming step further 
includes transforming second data from a second acoustic signal into the frequency domain, and 
wherein said frequency semblance value is additionally calculated from said second data 
10 8. The method of claim 5, further comprising: 

(c) graphically displaying a set of said frequency semblance values. 

9. The method of claim 5, further comprising: 

(d) generating a set of said frequency semblance values. 

1 0. The method of claim 9, further comprising; 

15 (e) applying a frequency filter to said data, said frequency filter corresponding 

to one or more ranges of frequencies, and wherein said frequency filter eliminates portions 
of said data that correspond to said ranges of frequencies. 

11. The method of claim 10, wherein said one or more ranges of frequencies are 
selected based upon predetermined patterns in said set of frequency semblance values. 

20 1 2. The method of claim 5, wherein a set of frequency semblance values are calculated, 

said set of frequency semblance values including at least one maximum semblance value, said 
maximum semblance value corresponding to a first peak frequency, a first peak slowness, and an 
energy value. 

13. The method of claim 12, wherein said at least one maximum semblance value 
25 corresponds to a set of high frequency semblance values, said set of high frequency semblance 

values corresponding to a characteristic of said first acoustic signal. 

14. The frequency domain method of claim 5, wherein said first acoustic signal is 
received at a plurality of depths in a borehole, said transforming step further including the 
transformation of data from each of said plurality of depths in said borehole for calculation of said 

30 semblance value. 

15. The frequency domain method of claim 14, wherein said set of frequency 
semblance values are calculated downhole in a borehole. 

16. The frequency domain method of claim 14, wherein said set of frequency 
semblance values are calculated at an above surface location. 
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17. The frequency domain method of claim 14, wherein a set of frequency semblance 
values are calculated and said set of frequency semblance values correspond to a first peak 
frequency, a first peak slowness, a peak energy value, and a peak shape. 

18. The frequency domain method of claim 17, wherein said first peak frequency, said 
5 first peak slowness, and said peak energy value are calculated by a downhole processor in a 

wellbore. 

19. The frequency domain method of claim 18, wherein at least one of said first peak 
frequency, first peak slowness, peak energy value, and peak shape are calculated at a sub-surface 
location to yield a calculated value, said at least one calculated value being stored downhole, said 

10 method further comprising calculating at least a second calculated value at an above ground 
location. 

20. A method for data compression, comprising: 

(a) recording a first waveform at a plurality of data points; 

(b) compressing said data points by a first protocol to derive a first compressed 

IS data; 

(c) compressing said data points by a second protocol to derive second 
compressed data; 

(d) comparing said first compressed data to said second compressed data to 
determine which would require less memory for storage. 

20 21. The method of claim 1 8, where said first protocol measures the difference between 

data points. 

22. The method of claim 19, wherein said first protocol measures the difference 
between measured differences. 

23. The method of claim 22, wherein said steps of compressing include compressing 
25 said data points by a third protocol to derive third compressed data, said first protocol being a bit 

encoding protocol, said second protocol being a difference encoding protocol, and said third 
protocol being a difference of the differences protocol. 

24. A frequency domain semblances method comprising: 

(a) generating a first acoustic signal at a first depth in a borehole; 
30 (b) receiving said first acoustic signal at a first set of acoustic receivers 

proximate said first depth to generate a first series of received acoustic waveforms; 

(c) transforming said first series of received acoustic waveforms into the 
frequency domain; 

(d) calculating a first frequency semblance value from said data, said frequency 
35 semblance value corresponding to a first frequency and a first slowness for said data; 
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(e) generating a second acoustic signal at a second depth in said borehole; 

(f) receiving said second acoustic signal at a second set of acoustic receivers 
proximate said second depth to generate a second series of received acoustic waveforms; 

(g) transforming said second series of received acoustic waveforms into the 
5 frequency domain; 

(h) calculating a second frequency semblance value from said data, said second 
frequency semblance value corresponding to a second frequency and a second slowness for 
said data; 

(i) correlating said first slowness and said second slowness with said first depth 
10 and said second depth. 

25. The method of claim 24, wherein said correlating step yields a slowness-depth log 
of said borehole. 

26. A method to combine time domain semblances, comprising: 

(a) calculating from a first plurality of waveforms a first set of time domain 
15 semblance values corresponding to a first depth in a borehole, each of said time domain 

semblance values in said first set corresponding to a time and a slowness; 

(b) calculating from a second plurality of waveforms a second set of time 
domain semblance values corresponding to said first depth in a borehole, each of said time 
domain semblance values in said second set corresponding to a time and a slowness; 

20 (c) combining a first time domain semblance value from said first set and a 

second time domain semblance value from said second set, said first time domain 
semblance value and said second time domain semblance value corresponding to the same 
time and slowness. 

27. The method of claim 26, wherein said combining step is multiplication. 

25 
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